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SUMMARY 

The paths of cloud droplets into two engine inlets have teen calcu- 
lated for a wide range of meteorological and flight conditions. The 
amount of water in droplet form ingested by the inlets and the amount 
and distribution of water impinging on the inlet walls are obtained from 
these droplet-trajectory calculations. In both types of inlet, a pro- 
late ellipsoid of revolution (10-percent thick) represents either part 
or al 1 of the forebody at the center of an annular inlet to an engine. 
The configurations can also represent a fuselage of an airplane with 
side ram-scoop inlets. The studies were made at an angle of attack of 
0°. The principal difference between the two inlets studied is that the 
inlet air velocity of one is 0.7 that of the other. The studies of the 
two velocity ratios lead to some important general concepts on water in- 
gestion in inlets. 


INTRODUCTION 

As part of a comprehensive program on cloud-droplet-trajectory com- 
putations, the NACA Lewis laboratory has calculated the paths of drop- 
lets ahead of and surrounding prolate ellipsoids of revolution moving 
through clouds (refs. 1 to 4) . The present report extends the trajec- 
tory computations with respect to prolate ellipsoids in an incompressible 
flow field to include the study of water ingestion in two engine inlets 
at 0° angle of attack. In both inlets, an ellipsoid represents either 
part or al 1 of the forebody at the center of an annular inlet to an en- 
gine, or some other air-processing device. In either inlet, the ellips- 
oid can also be assumed to approximate the fuselage of an airplane such 
that the inlets become side ram-scoop inlets. 

The configurations, which are described in detail in the following 
section, are geometrically simple in order to permit the studies de- 
scribed herein. Jn spite of the simplicity, the configurations are 



2 


NACA TN 3593 


reasonable approximations to those found on aircraft. The results of 
the studies on the two simple shapes lead to important general con- 
cepts on water ingestion into inlets of this type. 


DESCRIPTION OF CONFIGURATIONS 

In order to make a study of impingement on inlets, it is necessary 
to devise a model configuration for which the air-flow field can he oh- in 

tained and which reasonably represents the actual case of an inlet on o 

an airplane. The models chosen for this study are basically dependent 
on the air velocity flow field surround in g a prolate ellipsoid of revo- 
lution. Both inlet configurations consist of an annulus surrounding an 
ellipsoid of revolution as the forebody. 

The inlets differ somewhat in physical appearance as well as in 
the entrance velocity of the air. In one inlet configuration, the cowl, 
or outer wall, is infinitely thin and conforms to the shape of an air 
streamline. With this inlet there are many possible locations of the 
outer wall with respect to the ellipsoid forebody. A possible configura- 
tion is shown in figure 1. The cowl of the other inlet configuration 
has a slight thickness, as shown in figure 2. This outer wall has a 
fixed configuration with respect to the ellipsoid forebody AB. The 
In ner surface of the cowl and the inlet inner wall are straight and 
parallel. Either of the inlets shown in figures 1 and 2 may be a com- 
plete annulus surrounding the forebody or a sector of an annulus such 
as would be found in side ram-scoop inlets. 

The ellipsoid chosen for this study is 10-percent thick (fineness 
ratio of 10) . The problem is limited to an angle of attack of 0° be- 
tween the major' axis of the ellipsoid and the free-stream air, because 
the facile methods of calculating the flow field and trajectories are 
limited to 0°. 

Each of the two configurations has an associated inlet velocity 
ratio, which is defined as the ratio of the velocity at the inlet en- 
trance to the free-stream velocity. For configuration 1 (fig. l), the 
inlet velocity ratio is nominally 1.0. (Precise calculations show that 
the inlet velocity ratio may vary from 0.98 to 1.02, depending on the ■ 
location of the inlet entrance.) The inlet velocity ratio for config- 
uration 2 (fig. 2) is nominally 0.7 (0.714, actually). A particular 
air-flow field is associated with each velocity ratio. The methods used 
in calculating the two flow-field patterns produced a cowl and inner 
wall with a slight difference in geometrical shape for the two inlets. 

In spite of this, the difference in velocity ratio is not caused princi- 
pally by the difference in geometry, but rather by differences in flow 
requirements . As will be discussed further, the slight cowl thickness 
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shown in figure 2 will not appreciably affect the amount of water in- 
gested, thus permitting the comparison of the two inlets on the basis 
of velocity ratio. 


Inlet with Velocity Ratio of 1.0 

The droplet impingement on the cowl and the amount of water in- 
gested can be obt ain ed by extending the analysis of the trajectory cal- 
culations made for the results reported in reference 4. Droplet tra- 
jectories in space surrounding the 10-percent-thick ellipsoid were 
calculated in order to obtain the results reported in reference 4 on 
spatial variation of local liquid-water concentration. These same tra- 
jectories are not altered by an inlet such as shown in figure 1 if the 
air str eamlin es (or air-flow field) ahead of and surrounding the ellips- 
oid are not altered by the inlet. The placement of an inlet around the 
ellipsoid (or around a sector between two meridian planes) will not alter 
the air -flow field if the cowl is of negligible thickness, the cowl 
contour coincides with that of the air streamline on which it is super- 
posed, and the values of local air velocity around the ellipsoid are 
not changed. The provision that the local velocity is not changed im- 
plies that the air flow remains ideal and without viscosity or boundary 
layers at the inlet walls. 

The example shown in figure 1 illustrating the geometrical config- 
uration discussed in the preceding paragraphs has the forebody ABC de- 
fined by the equation z^ + 100r^ = 0.25. This is the equation for a 
10-percent-thick ellipse located with zero coordinates at the midpoint. 
The outer wall DE of the inlet is given by the equation for the air 
streamlines : 

( \ + 1 \ X I -L/2 

Ml - 1/ -*2 n 1 , /I + e\ 1 + 

x • Mr — ; - pTTa 

|(%) 2 (X - nfKX 2 - 1) (1) 

Symbols used herein are defined in appendix A. 

Air streamlines are shown as dashed lines in figure 3.- A much 
simpler equation of the form 


t-iTd- - !) ; 


z^ + cr^ = g (2) 

applies very well over the region between z = -0.4 and z = 0. The 
ordinate scale r in figure 3 is expanded to 4 times the abscissa 
scale z in an effort to obtain precision. This scale distortion is 
used in several subsequent figures presenting trajectory data. 



4 


NACA TN 3593 


Trajectories' of droplets are shown in figure 4. Droplets shown 
starting at infinity at points j and k impinge on the forebody at 
points j s and k s , respectively, and are assumed to be removed from 

the possibility of entering the inlet. The trajectory with starting 
point t is tangent to the forebody at t s and continues in space 

around the elliptical forebody. Information on droplet impingement be- 
tween r Q = 0 and the tangent trajectory is presented in reference 3. 

Information on the shadow zone formed between the tangent trajectory 
and the ellipsoid surface (from t s back to the ellipsoid midsection) 

is presented in reference 4. 

Three possible entrance walls are shown in figure 4. The choice 
of either entrance wall IK or GH or any other wall such as DE (fig. 4) 
that coincides with an air streamline (fig. 3) depends on entrance-size 
requirements. If an entrance wall such as <JK (fig. 4) meets the 
entrance-size requirements,, then all the droplets in the space between 
the tangent trajectory and the entrance lip trajectory l-l-p will be 

ingested. (This discussion applies only to clouds made of droplets of 
uniform size.) The space between the trajectories is bounded by the 
surfaces formed by rotating the trajectories in the r,z -plane about the 
axis of the ellipsoid. With the entrance wall placed at GH, trajec- 
tories such as n-n with starting ordinates r Q larger than those for 

m-mp will miss the forebody and entrance altogether. The liquid water 
in droplet form in the annulus space between trajectories m-mp and 
q-q will impinge on the inner surface of wall GH between z = -0.25 
and 0. The detailed methods for finding both the amount of water in- 
gested and the water distribution are shown with examples in a subse- 
quent section. 


Inlet with Velocity Ratio of 0.7 

The solution of the air-flow field for an inlet velocity ratio 
other than 1.0 was obtained in a median plane by adding a two-d i mensional 
perturbation field, representing the effect of the inlet, to the r- and 
z-camponents of the air velocity field surrounding a 10-percent- thick 
prolate ellipsoid of revolution. The addition of a two-dimensional 
perturbation field to a three- dim ensional field confined to one meridian 
plane is not strictly correct. This method of solution is justified 
only if the two-dimensional perturbation is not large and if the inlet 
is located where the radial component of flow is small, that is, not 
near the nose of the ellipsoid. The two requirements are complied with 
in the solution for this problem, because the strength of the perturba- 
tion reduces the velocity ratio to only 0.7 and the inlet is located 
outside the ellipsoid in a region of largest ellipsoidal diameters. 
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The perturbation field used in this study is derived and discussed' in 
appendix B. The boundary conditions at the surfaces of the forebody 
and inlet walls are satisfied, because these surfaces represent air 
streamlines after the perturbation velocities are added to the basic 
flow-field velocities of the 10-percent- thick ellipsoid. The result ing 
configuration shown in figure 2 has an attendant flow field that was 
used in the calculation of the trajectories. "Unlike the preceding case 
for a velocity ratio of 1.0, the cowl is fixed in its location with re- 
spect to the inner wall and centerline of the forebody, and only one 
configuration is associated with this inlet velocity ratio of 0.7. 

The forebody and inlet inner wall ABC (fig. 5) are defined para- 
metrically by 


-z = 0.00943 (<P + 0.9956^) + 0.2780 

r =? 0.00943 (0.1 + 0.09986^) + 0.0418 

for values of <P between 0.6 and -29.5. This curve approaches 
r = 0.041 between z = -0.25 and z = 0, and very closely approximates 
an ellipse between z = -0.30 and z = -0.47. In order to simplify the 
concepts and the presentation of the data, the forebody AB in figure 5 
and all subsequent figures illustrating this inlet is the 10-percent- 
thick ellipse defined by the equation z 2 + 100r 2 = 0.25. The outer 
wall DE of the inlet is defined by a section of the ellipse 

z 2 + 133r 2 = 0.75. The inner surface DF of the cowl is the line 
r = 0.0714. 

As in the case for a velocity ratio of 1.0, the amount of water in- 
gested into the inlet is the water contained in the space between the 
tangent trajectory (t-t s , fig. 5) and the trajectory that intersects the 

entrance lip on the cowl of the inlet (trajectory Z-7 p , fig. 5). 


DETERMINATION OF DROPLET PATH 

The method for the determination of the droplet path Is essen- 
tiality the same as is described in reference 1. A differential analy zer 
of the mec hani cal- anal og type was used to solve the equations of motion 
of the droplets in the airstream for various values of the parameter 
l/K between 0.1 and 90. The inertia parameter K is a measure of the 
droplet size, the flight speed, the size of the basic ellipsoid, n.nd 
the viscosity of the air, in the form 


K = 1.704X10 


-12 d% 
|iL 


( 3 ) 


The density of water, which is expressed as part of the conversion 
factor, is 1.94 slugs per cubic foot. 
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For each, value of the parameter l/K, a series of trajectories was 
computed for each of several values of free-stream Reynolds number Reg . 

The free-stream Reynolds number with respect to the droplet diameter is 
defined as 


<lP a u 

Re 0 = 4. 813X10" b — — 




In order to make these dimensionless parameters more physically signifi- 
cant in the following discussion, seme typical combinations of K and 
Reo are presented in table I for various lengths and velocities of the 
basic e lli psoid, droplet sizes, and flight pressure altitudes and tem- 
peratures. A procedure for rapid calculation of l/K and Reg from 
practical flight conditions is given in appendix B of reference 5. 

Before the integration of the equations of motion (ref. l) to ob- 
tain the trajectories could be performed with the differential analyzer, 
the initial velocity of the droplets had to be determined at the point 
selected as the starting position. In addition, since the starting po- 
sition must be selected at a finite distance ahead of the ellipsoid, it 
was necessary to make a correction to this starting ordinate in order 
to obtain the corresponding starting ordinate rg at z = -<*>. The 

procedure for obtaining the starting conditions for the analog is the 
same as that discussed in detail in reference 1. 


RESULTS AND DISCUSSION 
Inlet with Velocity Ratio of 1.0 

Several sets of trajectories for different combinations of Reg 

and l/K are shown in figure 6. Same of the many possible wall loca- 
tions are shown in figure 6 as dashed lines, which are duplications of 
the air streamlines shown in figure 3. The cowl lip may be located 
anywhere along the air streamline. The trajectories are shown as 
medium-weight solid lines. The order of presentation of figures 6(a) 
to (e) is in accordance with the amount that droplet trajectories dif- 
fer from air streamlines. That is, in figure 6(a) the difference be- 
tween a droplet trajectory and a superposed air streamline (particu- 
larly regarding the slope of each at a given point in space) is greater 
than in figure 6(e). In figure 6(e) the directions of the droplet tra- 
jectories are everywhere nearly coincident with air streamlines. This 
order of presentation also shows the correlation between droplet im- 
pingement on the cowl and the parameters Reg and l/K. The impinge- 
ment in figure 6(d) is less than in figures 6(a) or (b) . 
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The significance in the variation of l/K and ReQ is seen in 

table I. Some of the approximate possible flight and atmospheric com- 
binations that could be applied to figure 6 are given in the following 
table: 


Flight 

speed, 

mph 

Pressure 
altitude , 
ft 

Droplet 

diameter, 

<1, 

microns 

Ellipsoid 
major 
axis, L, 
ft 

Figure 

300 

5,000 

1000 

50 

6(a) 


15,000 

80 

29 

6(b) 



20 

10 

6(c) 




27 

6(d) 

50 

25,000 

10 

3 

6(e) 


Droplet ingestion . - The amount of water ingested by an inlet is 
the water in droplet form through an annular area of space (perpendicu- 
lar to the major axis of the ellipsoid) between the inlet lip and the 
body of the ellipsoid. For example, if the inlet cowl were located at 
JK in figure 4, water would be ingested through an annular area of 
space of width BJ. For the case of a velocity ratio of 1.0, wherein 
any number of locations of the cowl are possible, the amount of water 
ingested must be presented as a function of inlet entrance location as 
well as Reg and l/K. 

If the point of tangency t g of the tangent trajectory occurs 

outside the entrance, such as in figure 4 (i.e., the z-position of the 
point of tangency is upstream from the z-position of the entrance lip), 
the amount of water ingested is obtained frcm the following equation: 

W-l = 0 . 33jtwoL^U(rQ^ ^ - r^ taQ ) (5a) 

where is the free-stream radial distance at z = -<*> of the tra- 

jectory that touches the entrance lip, and r o_,tan ^- s free-stream 

radial distance of the tangent trajectory. If the tangency point is 
located downstream of the entrance lip (i.e., inside the entrance, see 
fig. 6(a)), the water ingested is obtained from 

% = 0.33*w 0 L 2 (r^j - r^ B ) (5b) 

where tq b is the radial distance at z - -<*> of the trajectory that 

impinges on the el li psoid at the same z-position as the inlet opening, 
such as point B in figure 4. The z-position of the tangent trajectory 
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can "be obtained from figure 8 of reference 4. The amount of water in- 
gested for various radial locations of the lip can be found frcrn the 
data presented in figure 4 of reference 4 if the inlet lip is located 
at z=0, -0.25, or -0.5. Intermediate locations of the entrance lip 
may be treated by cross-plotting the data in reference 4. 


The following example illustrates the use of figure 4 in refer- 
ence 4. Assume flight and atmospheric conditions such that Beo = 128 
and l/K = 45. These conditions correspond to an approximate flight 
speed of 300 miles per hour, cloud-droplet diameter of 20 microns, 
ellipsoid length of 27 feet, and a pressure altitude of 15,000 feet. 
Further, assume that design criteria governing the inlet spacing area 
place the inlet lip at r = 0.08 and z = -0.25. Figure 4(h) in ref- 
erence 4, which is reproduced herein for convenience as figure 7, is 


used to obtain the required values of r; 


0,1 


and 


'0,tan' 


From either figure 7 herein or figure 4(h) in reference 4, the 

value of r 2 = 0.00477 is found from the curve for l/K = 45 for an 
0,7 1 

abscissa reading of r^ = (0.08)^ = 0.0064. The value of r^ can 

be obtained from the same figure by reading the bottom end of the same 
curve. For this example, r 2 = 0.000045. The only remaining un- 
known for the solution of equation (5a) is the liquid- water content w Q , 
which can be obtained from meteorological data. The bottom ends of the 
curves in figure 4 of reference 4 automatically give either the value 
of r 2 i _ an or r 2 B , whichever applies. If the bottom ends of the 

curves touch the dashed line labeled "ellipsoid surface," Tq -g is ob- 
tainedj whereas rQ,tan is obtained from the curve ends that do not 
touch the dashed line. 


If the cowl lip is located at z = -0.5 and the cowl coincides 
with an air streamline, the water ingested can be obtained from the data 
in figure 4(m) of reference 4. The procedure for determining the water 
ingested is the same as in the preceding example, except that 

W ± = 0 . 33nwL^UrQ^ ^ (5c) 

because the lower limiting trajectory coincides with the z-axis. 

The method described for obtaining the amount of water ingested 
pertains in principle to a simple configuration in which a portion of 
an air streamline is replaced by a wall of negligible thickness. In an 
actual model, the wall has thickness and probably deviates from an air 
streamline. The method still applies if the wall thickness is not 
large compared with the ellipsoid thickness, which is usually the case, 
and the centerline of the wall shape approximates an air streamline. 


in 
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In other words , the results obtained by this method are not sensitive 
to small deviations from the ideal model presumed in these studies, be- 
cause small amounts of wall thickness and wall deviation from the origi- 
nal shape of the air streamline introduce only small changes in the air- 
streamline pattern. Furthermore, because of the droplet inertia, the 
change in droplet trajectory is even less than the change in air- 
streamline pattern. 

Total impingement on inlet walls . - The amount of water impinging 
on the inlet walls is a function of the inlet lip position as well as 
of the values of l/K and ReQ. If the point of tangency t s is out- 
side (upstream from) the entrance opening (fig. 4), there is no impinge- 
ment on the inner wall of the inlet. The z -position of the point of 
tangency t s , which corresponds to zero value of shadow zone, is given 

in figure 8 of reference 4. 

If the point of tangency occurs inside the entrance, the amount of 
water impinging on the inlet inner wall (i.e., on the ellipsoid surface 
downstream of the inlet opening) can be found from the data in refer- 
ence 3. The water impinging on the ellipsoid surface between the en- 
trance opening and the point of tangency t s on the inlet inner wall is 

W B-t s = 1 * 04w 0 L2 ( r 0,tan " r 0,B) ( 6 ) 

where 1.04 is the combination of constants in equation (5) . The value 
of Tq tnr| Is available directly in figure 5 of reference 3. The value 

for r 0 p can be obtained from figure 4 of reference 3. A correlation 

between z and the distance S along the ellipsoid surface for use in 
connection with figure 4 of reference 3 is available in figure 3 of the 
same reference. 

Whether impingement occurs on the outside surface of the cowl de- 
pends primarily on the values of l/K and ReQ. For example, the tra- 
jectories shown in figure 6(a) are nearly straight lines. For the con- 
ditions represented by figure 6(a), impingement occurs only on the 
outside surface of the cowl. For conditions represented in figure 6(b), 
the amount of impingement on the outside and inside cowl surfaces de- 
pends on the r- and z-position of the inlet lip. For conditions repre- 
sented in figures 6(c) to (e), impingement occurs only on the inside 
surface. 

Figures 6(a) and (e) represent the two extremes of trajectory 
shape. Combinations of values of l/K less than l/30 with ReQ less 

than 8192 result in trajectories at least as straight as those shown in 
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figure 6(a) . Combinations of values of l/K greater than 114 with 
ReQ greater than 8 result in trajectories that follow the air stream- 
lines at least as well as those shown in figure 6(e) . For values of 
l/K less than 3, the impingement is principally on the outside surface 
of the cowl between z = -0.5 and -0.25. For values of l/K greater 
than 3, the impingement is principally on the inside surface spread 
thinly over the entire surface. 

Impingement distribution on inlet walls . - No attempt is made to 
give quantitative values of the amount of water impinging on the cowl 
nor quantitative values of the distribution of impingement. Both the 
amount of impingement and the distribution are very sensitive to small 
changes in the shape of the cowl, particularly for conditions in which 
the droplet trajectories are nearly like air stre amli nes (figs. 6(c) to 
(e) ) . The use of figure 6 as a guide will lead to useful approximate 
results when the cowl has some thickness and the centerline can nearly 
be superposed on an air streamline. For flight and meteorological con- 
ditions represented by figures 6(c) to (e), cowl thickness will bring 
about impingement on the entrance lip. The collection on the lip will 
be very high. For these same flight and meteorological conditions, im- 
pingement on the outside surface is hardly perceptible. On the inside 
surface the impingement is slight but spreads thinly for a considerable 
distance downstream of the lip. 

Impingement for conditions shown in figure 6(a) can be obtained 
very easily from the geometry, because the trajectories are nearly 
straight lines. The most difficult problem is presented by conditions 
shorn in figure 6(b) . It is helpful to know that for the conditions 
shown in figure 6(b) a cowl with thickness does not disturb the tra- 
jectories, because the inertial forces of the droplets are relatively 
large compared with the forces in the airstream available to change 
their paths. Therefore, inlet walls with small amounts of thickness 
can be superposed directly on figure 6(b), and the impingement and dis- 
tribution can be obtained from the resulting geometry. 

Variation of liquid water across inlet opening . - The installation 
of small protruding objects that do not alter the air flow appreciably 
is often required near the inlet opening. Occasionally, coarse screens 
are also used, A knowledge of the rate of ice formation on these ob- 
jects is very valuable. As with the ellipsoids reported in references 
2 and 4, a radial variation in the local liquid-water concentration 
exists in the Immediate vicinity of the elliptical forebody. Thus, the 
local concentration of water varies radially at the inlet opening. The 
procedure for determining the spatial variation of liquid water is de- 
scribed in detail in references 2 and 4. 

20-Percent-thick el 1 i psoid . - The preceding results and discussion 
have been presented for a 10-percent-thick ellipsoid. The same 
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discussion and methods apply to a 20-percent-thick e lli psoid. The cor- 
responding data are available in reference 2, which presents figures 
that correspond in notation and ordinates to those cited in reference 
4. The same correspondence of figures exists between references 3 and 
1, which are required in order to determine the impingement on the 
inner wall (i.e., the ellipsoid). A figure of the air stre amlin es 
surrounding a 20-percent-thick ellipsoid is given herein as figure 8. 
The air streamlines can also be calculated from equation (l) . Figure 
6 herein may be used as a guide for estimating the impingement and dis- 
tribution on the outer wall of the 20-percent-thick body as well as the 
10-percent- thick ellipsoid. The dual use of figure 6 is possible be- 
cause,, relative to the air streamlines , the droplet trajectories are 
very similar for the two bodies. 


Inlet with Velocity Ratio of 0.7 

The effect of changes in l/K and Reo on droplet trajectories 
is illustrated in figure 9. The same range of flight and meteorologi- 
cal conditions covered by figure 6 is given in figure 9. 

Droplet ingestion . - The amount of water in pounds per hour flow- 
ing in the region bounded by the tangent trajectory (t-t s , fig. 5) and 

the upper limiting trajectory (Z-7p, fig. 5) is given by the expression 

W ± = 0.33w 0 UAC t (7) 


The derivation of this equation is based on the law of conservation of 
matter and is similar to the derivation of equation (4) in reference 2. 

The total concentration factor C^. may be considered an ingestion 

efficiency (scooping ratio), since it is the ratio of the area at 
z = - ® between the tangent and upper limiting trajectories to the area 
at the inlet opening. The term C-j. is defined as 

2 2 

G _ 0,Z ~ 0,tan 

t (0.0714) 2 - (0.04) 2 


where 


r 0 ,1 
tory, and r. 


is the ordinate at z = -«° of the upper limiting trajec- 

L 0 tan or( 3i nfl ‘* :e a i infinity of the tangent trajectory. 

The values in the denominator are the values of r of the inner wall and 
cowl at the opening (points B and D, fig. 2). This definition of the 
total concentration assumes that the area of interest A, in square 
feet, is a sector of an annulus with the area perpendicular to the major 
axis of the ellipsoid. The values of are given in figure 10. 



12 


NACA TN 3593 


The values tq j of the upper limiting trajectories and r o ; tan 

of the tangent trajectories are given in figure 11. The trajectories 
near the body are deviated more from a straight line than those near 
the lip of the entrance (see fig. 9), so that the ordinates at z = 
of the upper limiting trajectories are not affected by changes in ReQ 

and l/K as much as the tangent trajectories. 

Variation of liquid water across inlet opening . - The local weight 
flux of water in droplet form is derived in reference 2 as 

F = 0.33wqUC, lb/(hr) (sq. ft) (9) 


where 



The element of area implied by the equation is a unit perpendicular to 
the major axis located in the plane of the inlet opening. The local 
concentration factor C at a point r in space at z = -0.30 and -0.25 
is obtained from figure 12. The largest values of concentration factor 
are near the forebody, just above the shadow zone. The variation in 
local concentration radially across the inlet is very similar to the 
spatial variation for the inlet with a velocity ratio of 1.0. These 
variations are described in references 2 and 4. 

Shadow zone . - A region of zero concentration exists adjacent to 
the surface of the inner wall. The thickness (r-r s ) of this shadow 

zone at each z -position on the body for various values of ReQ and 

l/K is given in figure 13. The magnitude of the shadow zone depends 
on both Re 0 and l/K, as is the case with the ellipsoids described in 

references 2 and 4. For the inlet with a velocity ratio of 0.7, the 
thickness of the shadow zone becomes appreciable inside the inlet, be- 
cause the inner wall becomes a straight line parallel to the major axis. 

Total impingement on cowl . - Although the wall thickness does not 
appreciably affect the total water ingested, the wall shape is impor- 
tant in determing impingement on the wall. With the parti cular value 
of inlet velocity ratio studied (0.7), no impingement occurs on the 
outer surface of the cowl for values of l/K larger than 1. This 
means that for most flight and atmospheric conditions (see table I) 
only large drops will impinge on the outer surface of the cowl. 
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The amount of water in droplet form impinging on the inner surface 
of the cowl is 

W c = 0.333W 0 UL 2 (rg^_ lp - r^_ q ) (10) 

O r\ 

The values of ( r o^Z-2p ” r 0,q-q) required, in equation (10) are given in 

figure 14. The impingement study is confined, to the surface between 
points D and F, because trajectories were not calculated beyond the mid- 
section of the ellipsoid (z = 0) . 

Impingement distribution on cowl . - The rate of water impingement 
per unit area on the inner surface of the cowl between points D and F 
can be found by applying the values of f3 c given in figure 15 to 


W M = 2.3v 0 UP c (11) 

This equation is derived from the concept that the flux of water in 
droplet form through an annular area of space width r 2 - rp impinges 

on the outer wall surface between points S c ^ 2 and- S c 

rt(r 2 " r l ^0 ^0 

W P,c “ 0,33 W 0 U 2jt(0.0714)(S 2 - S x ) c “ 2 ' 3 W 0 U AS^ 


where , in the differential limit, 

dS c 


Pc 


The local impingement efficiencies P c are given in figure 15 as 
functions of the surface distance S c , starting from the inlet lip at 
point D. For all values of 1/k and ReQ the maximum impingement 

rate occurs at the entrance lip. Although the high rate occurs between 
the entrance lip and approximately S c = 0.02 (ratio of actual distance 

to major axis), the remainder of the surface is exposed to impingement. 
The impingement distribution on the inner surface of the cowl is ob- 
tained directly from figure 15, because Wp^ c is directly proportional 

to p c . The integral of P c between the cowl lip and z = 0 is equal 
to the total impingement on the cowl given in the preceding section. 


COMPARISON OF ENTRANCES WITH VELOCITY RATIOS 1.0 AND 0.7 

A comparison of the two types of entrances discussed herein can be 
done in a limited manner only. The configuration with the inlet 
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velocity ratio of 1.0 must "be confined, to that configuration which most 
nearly corresponds to the fixed, geometric shape of the 0.7-velocity- 
ratio inlet. An air streamline (line DE, fig. 4) exists for the inlet 
with a velocity ratio of 1.0 that superposes very well on the w all DE 
(fig. 5) of the 0.7-velocity- ratio inlet. An inlet composed of ellipti- 
cal forehody ABC and outer wall DE, shown in figure 4, is compared with 
an inlet configuration shown in figure 5 . The principal geometric dif - 
ference in the two configurations is that the outer wall DE in figure 4 
has negligible thickness; whereas, the outer wall in figure 5 has thick- 
ness EP at z = 0. In spite of the thickness, EF the outer wall of 
the 0.7-velocity- ratio inlet is very thin, approximating a thin wedge 
(in sectional view) with a sharp edge at the entrance point D. The 
relative geometric proportions of the 0.7-velocity-ratio inlet are best 
illustrated in figure 2. (As was previously stated, the ordinate scale 
in figure 5 is expanded with respect to the abscissa scale.) The thick- 
ness of the wall is not enough to be of consequence in the influencing 
of the air-flow field ahead of the inlet. The inlet inner walls are 
also shaped somewhat differently; the inlet opening areas are the same 
and in both inlets remain approximately constant from the opening 
through to z = 0. In spite of .the smal 1 differences mentioned, the 
two inlet configurations are enough alike to permit the comparison of 
total water ingested and outer- wall distribution for the two velocity 
ratios . 


Total Water Ingested 

In order to facilitate the presentation of results, the comparison 
is made at several specific flight and atmospheric conditions. The 
conditions chosen Eire the same as those represented by the trajectories 
of figures 6 and 9. These conditions are representative of a wide 
range of values of l/K and Reg. However, a danger must be noted 

arising from the choice of only a few sets of conditions for comparison: 
the relative importance of the different factors, such as flight speed, 
droplet diameter, and so forth, being compared may change when other 
numerical values are chosen. 

The comparative results are presented in the following tabulation 
as ingestion efficiencies. The ingestion efficiency (or scooping ratio) 
is the ratio of the water ingested to the total water in an annular 
space which has a cross-sectional area equal to that of the inlet and a 
volume in free space equal to that swept out by the inlet: 
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Flight 

speed, 

mph 

Pres- 

sure 

alti- 

tude, 

ft 

Droplet 

diam- 

eter, 

<1, 

microns 

Ellips- 

oid 

major 

axis, 

L, 

ft 

Recip- 
rocal 
of in- 
ertia 
param- 
eter, 
l/K 

Free- 

stream 

Reynolds 

number, 

Re 0 

Ingestion 

efficiency, 

percent 

1.0- 

Velocity- 

ratio 

inlet 

0.7- 

Velocity- 

ratio 

inlet 

300 

5,000 

1000 

50 

l/30 

8192 

100 

100 


15,000 

80 

29 

3 

512 

102 

94 



20 

10 

15 

128 

99 

89 




27 

45 

128 

97 

83 

50 

25,000 

10 

3 

114 

8 

94 

78 


The ingestion efficiency includes the water that impinges on the 
inner wall. For l/K = l/30, ReQ = 819 2 , and a 1 . 0- velocity-ratio in- 
let, 8 percent of the water passing through the opening strikes the in- 
ner wall (impinges on ellipsoid surface beyond the opening) . For a ll 
other conditions tabulated, all or almost all the water ingested does 
not strike the inside surface of the cowl up to the midsection (z = 0) . 
For l/K = 3, Reg = 512, and a 1.0- velocity- ratio inlet, a slight con- 
centration effect is present, which results in an ingestion efficiency 
of 102 percent. The ingestion efficiency of the inlet with a velocity 
ratio of 1.0 is over 90 percent and less than 104 percent for a wide 
range of conditions listed in table I. 

The ingestion efficiency of the inlet with a velocity ratio of 0.7 
is always over 75 percent and never over 100 percent for the conditions 
listed in table I. The ingestion efficiency for this inlet is found 
directly from figure 10, which presents the total concentration factor. 
Ingestion efficiency is 100 times larger than the dj. factor of figure 

10. For comparable flight and atmospheric conditions, somewhat less 
water Is ingested at 0.7 velocity ratio than at 1.0, except for rain 
drops. 


Distribution of Water on Inner Surface of Cowl 

A detailed quantitative comparison of the water distribution has 
little practical application, because small changes in shape of the 
walls affect the distribution appreciably. A reasonable qualitative 
comparison can be made easily by examining figures 6 and 9; for ex- 
ample, for 1 /K = l/30 and ReQ = 8192 (figs. 6(a) and 9(a)), the tra- 
jectories for both inlets are practically straight lines, with consequent 
slight impingement on the outer surface and no impingement on the inner 
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surface of the cowl. For other combinations of l/K and ReQ in fig- 
ures 6(h) to (e) and 9(h) to (e) , in which there is some Impingement on 
the inner surface of the cowl for hoth inlets , the distribution of this 
impingement differs between the two inlets. For the inlet with a veloc- 
ity ratio of 1.0, the trajectories are nearly parallel to the cowl (lo- 
cated at DE, fig. 4) . The smalJ amount of impingement is spread fairly 
evenly from points D to E of the wall surface. For the inlet with a 
velocity ratio of 0.7, the distribution is much higher near the inlet 
opening (see fig. 15 for values of local impingement efficiencies) . 

Thus, lowering inlet velocity ratio from 1.0 to 0.7 tends to distribute 
the impingement less evenly on the inner surface of the cowl with greater 
concentration near the opening. 

For both inlets in an actual cloud in which a distribution of drop- 
let sizes is present, the area some distance beyond the lip (say 
S c >0.1) can be expected to build up unevenly with ice if this area is 

unprotected. The reason for this expected jaggedness in the ice build- 
up is that the droplets strike the wall surface at a small angle (tra- 
jectories nearly parallel to surface, see figs. 6 and 9). The points of 
impingement of the small droplets in the cloud are influenced by small 
irregularities of the wall surface and by small-scale turbulence near 
* the wall. The start of ice build-up at a point of wall irregularity 
produces an unstable condition that promotes further growth of ice on 
that particular point and produces further turbulence for other growths. 


Shadow Zone 

The shadow zone is increased considerably in the 0.7-velocity- 
ratio inlet as compared with the 1.0- velocity- ratio inlet, principally 
because of the change in physical geometry of the inner wall rather than 
the change in velocity ratio. The shapes of the tangent trajectories, 
which determine one boundary of the shadow zone, are very similar in 
the two inlets for comparable flight and atmospheric conditions. 

In both inlets and for some flight and atmospheric conditions the 
thickness of the shadow zone becomes appreciable inside the inlet. If 
possible, consideration of the shadow zone should be given in the posi- 
tioning of boundary- layer scoops. The thickness presented in figure 13 
herein, in figure 8 of reference 4, and in figure 9 of reference 2, is 
for 0° angle of attack. A boundary- layer scoop designed according to 
figure 13 to be in the shadow zone might be subjected to heavy concen- 
tration when the aircraft is yawed or pitched, because the thickness of 
the shadow zone is very sensitive to changes in angle of attack. 
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CONCLUDING EEMAEKS 

The data presented herein can be-used as a guide in the evaluation 
of water ingestion in inlets. The studies on the two simple shapes at 
0° angle of attack lead to some general concepts on water ingestion in 
inlets. The amount of water ingested is not sensitive to sma l 1 changes 
in shape of the outer wall. The placement of a wall in the air flow 
will usually disturb the air flow ihore than the droplet trajectories. 

The impingement on the cowl (i.e., amount and distribution) is quite 
sensitive to the physical shape and surface condition of the wall. In 
most applications impingement on the interior surface of the cowl can 
be expected for a considerable distance inside the entrance, because 
the forebody deviates the trajectories away from the body and toward 
the cowl. 

The use of screens and boundary- layer-removal scoops at the en- 
trance requires careful design because of the existence of the shadow 
zone and regions of high concentration. The variation in concentration 
of water across an inlet shows that a screen can build up ice unevenly 
in such manner as to result in large radial -flow distortions. 

Although only one inlet configuration was studied at an inlet ve- 
locity ratio of 0.7, the general concept that lowering the inlet veloc- 
ity ratio lowers the ingestion efficiency can be applied to a large 
variety of configurations. This concept is evident from a study of the 
manner in which air streamlines change when the inlet air velocity ratio 
is decreased for a particular configuration. Since most aircraft (par- 
ticularly interceptor types) should be designed with protection against 
ice formations during loitering, climb, and let-down, a design based on 
water- ingestion data for an inlet velocity ratio of 1.0 will be adequate 
for lower velocity ratios. A lowering of the velocity ratio will also 
result in higher impingement near the inlet opening on the inner surface 
of the outer w all . 

As can be seen from figures 6 and 9, no impingement occurs on the 
outer surface of the cowl over a large range of meteorological and 
flight conditions. High local impingement efficiencies can be expected 
on the cowl lip for al 1 entrance velocity ratios. 

During the major portion of most flight plans the inlet velocity 
ratio does not exceed 1.0 and does not decrease below 0.7. For most 
inlet designs, the accuracy of the methods presented herein for esti- 
mating amount and distribution of impingement is within the accuracy of 
knowledge of the meteorological factors involved, such as liquid-water 
content and droplet size. 

As was stated early in this report, droplets that impinge on the 
forebody are assumed to be removed from the possiblity of entering the 
inlets. In reality, the rater deposited on the forebody may ran back 
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over the "body surface and into the inlet when the ice is melted by ther- 
mal anti -icing devices. A discussion of runback is beyond the scope of 
this report. The conditions resulting in runback and the effect of run- 
back on screens, boundary- layer scoops, and other devices in the inlets 
should be carefully considered in inlet designs. 

Although the calculations were made for incompressible flow, they 
should be applicable throughout the subsonic region because of the small 
effect of compressibility on droplet trajectories (ref. 6) and the high 
flight critical Mach number of the configurations studied. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, November 2, 1955 
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APPENDIX A 
SYMBOLS 

The following symbols are used in this report: 

A area, sector of annulus in plane perpendicular to major axis, 

sq. ft 

C local concentration factor, dimensionless 

£). total concentration factor, dim ension l ess 

c,g constants for size of inlet (see eq.. (2)) 

d droplet diameter, microns 

F weight flux of water, lb/(hr)(sq. ft) 

K inertia parameter, 1. 704X10"^^ d^u/pL, dimensionless (eq.. (3)) 

L major axis of ellipse, ft 

Reg free-stream Reynolds number with respect to droplet, 

4.813X10 -6 dp u/p, dimensionless (eq.. (4)) 
el 

r,z cylindrical coordinates, ratio of actual length to major axis, 
dimensionless 

S distance along surface of ellipsoid measured from stagnation 

point, ratio of actual distance to major axis, dimensionless 

S c distance along inner surface of cowl measured from point of in- 

let opening, ratio of actual distance to major axis, 
dimensionless 

U free-stream airspeed, mph 

rate of water impingement on inner wall, Ib/hr 
s 

W c rate of water impingement on cowl, lb/hr 

rate of water ingestion, lb/hr 

Wo n local rate of water impingement on surface of cowl, 

p '° ib/(to)( sl ft) 
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w liquid-water content in cloud, g/cu m 

a l/4 focal distance of ellipsoid 

P c local impingement efficiency on cowl, dimensionless 

£ eccentricity of ellipse defined by a meridian section of ellips- 

oid of revolution 

X, prolate -elliptic coordinates 

p. viscosity of air, slugs/ (ft) (sec) 

p a density of air, slugs/cu ft 

potential function 
ijr streamline function 

Subscripts: 

B point of inlet entrance on ellipsoid surface 

l lip 

r radial component 

s shadow 

tan tangent 

z axial component 

0 free- stream conditions * 
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APPENDIX B 


DERIVATION OF PERTURBATION FIELD FOR 0 . 7- VELOCITY-RATIO INLET 

The air-flow field for the configuration shown in figure 2 with an 
inlet velocity ratio of 0.7 was obtained by adding a two-dimensional 
perturbation field to the basic flow field about a 10-percent-thick 
ellipsoid of revolution. The perturbation field had the configuration 
described in reference 7 and mentioned in reference 8 under the name 
"Borda mouthpiece." This perturbation field is expressed by the com- 
plex function 


= 0* + e 0 * 


(Bl) 

where 

• 


£* = z* + ir* 


(B2) 

0* = <P + i\|r 


(B3) 

and the asterisk denotes special coordinates 
Borda mouthpiece and differing from those of 

applicable only to the 
the basic ellipsoid. 


The rectangular coordinates are given parametrically by 


z* = <p + e'P cos t 


(B4) 

r* = t + e^ sin t 


(B5) 


The application of the Cauchy-Riemann conditions 



= 

Sr* 


* _ _ Sjf 

^ “ Sr* " " Sz* 


results in the velocity components : 



1 + e^ cos t 
1 +-2e^ cos ijr + e^ 


j*. e^ sin if 

cp 1 2T 

1 + 2e cos f + e 


(B6) 


(B7) 
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In order to give the Inlet proper position and size with respect to 
the basic ellipsoid, it is necessary to perform a translation and change 
of scale according to the following transformation: 

z = 0.00943z* + 0.2780 = 0.00943 (<P + e 9 cos t) + 0.2780 (B8) 

r = 0.00943r* + 0.0418 = 0.00943 (t + e^ sin t) + 0.0418 (B9) 

A z,r point described by these equations corresponds to the same 
point in space described by the z,r coordinates used for the basic 
ellipsoid. The transformation leaves the velocity components unchanged., 
but they must be expressed in terms of the free-stream velocity U and 
inlet velocity u as: 




1 + e^ cos 


u. 


r,p 


1 + 2e^ cos ijr + e 


e^ sin i|f 

1 + 2e^ cos Tjr + e^ 


ijr /. u\ 0.5 (l + e^ cos ilr) 

2(p \ ” U/ , _ cp , 2 <P 

_ 1 * ' 1 _L On ™ nnn iff _1_ n 




1 + 2e v cos i(f + e 

0.3e^ sin ijc 
1 + Ze* cos 


(BIO) 

(Bll) 


for an inlet velocity ratio of 0.7. The subscript p denotes the per- 
turbation to be added to the field about the ellipsoid. 

The transformation given in equations (B8) and (B9) approximately 
matches one streamline of the perturbation field to the surface of the 
basic ellipsoid. Adding the perturbation velocity components in equa- 
tions (B10) and (Bll) to the velocity components of the basic 10- 
percent -thick ellipsoid gives one streamline that is the forebody ABC 
shown in figure 2 and another streamline that is the inlet outer wall 
FOE. Also, when the perturbation velocities obtained from equations 
(B10) and (Bll) sire added to the basic components, the resulting veloc- 
ity at the center of the inlet opening is 0.714 with respect to the 
free-stream velocity U, but 0.7 with respect to the local basic- 
ellipsoid velocity. 

The velocity components of the basic flow field may be obtained 
from equations given in reference 1 and plotted in figure 2 of refer- 
ence 3. 
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TABLE I. - RELATION OP DIMENSIONLESS PARAMETERS TO ELLIPSOID SIZE AND ATMOSPHERIC AND PLIOHT CONDITIONS 
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ft 
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Temperature, 

op 





u, 
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ai- 


20 

1 
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Re 0 

K 

1/K 

Re 0 

/ 

K 

iA 

R«0 

K 

1/K 

Cloud 

50 

10 

3 

14,7 

0.008123 

123.1 

10.72 

0.008383 

119.3 

7.836 

0,008793 
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drop- 



30 
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50 

3 
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53.62 
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30 
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100 
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60 

176.4 

.01169 
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64.4 

0 . 001509 
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47.01 
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0.005846 
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300 

176.4 

.001949 

513.1 
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20 

10 
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10.26 
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50 

294.1 

.01949 

51.31 

214.5 

,02012 
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30.18 
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.06840 
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Figure 6. - Trajectories far inlet with velocity ratio of 1.0. 
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(b) Reciprocal of inertia parameter, 3; free-stream Reynolds number, 612. 

Figure 6. - Continued. Trajectories for inlet with velocity ratio of 1.0. w 
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(c) Reciprocal of inertia parameter,* 15; free-etream Reynolds number, 120. 
Figure 6. - Continued. Trajectories for inlet with velooity ratio of 1.0. 
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Axial dlatanoe, b 

(d) Reoiprooal of inertia parameter, 46 j free-strean Reynold® number, 1ES 
Figure 6, - Oontinued. TrajeotarieB for inlet with veloolt? ratio of 1.0 
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(e) Reoiprooal of inertia paronetar, 114; fraa-atreara Raynolda nunber, 6. 
Figure 6. - Concluded, Trajectories for inlet with velocity ratio of 1.0. 
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Axial distance, z 

(b) Reciprocal of inertia parameter, 3j free -stream Reynolds number, 512. 
Figure 9. - Continued, 'trajectories for inlet vitb velocity ratio of 0.7. 
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(d) Reciprocal of inertia parameter, 45j free-stream Reynolda number, 128. 

Figure 9. - Continued. Trajectories for inlet vith velocity ratio of 0.7. 
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(e) Reciprocal of inertia 
Figure 9. - Concluded. Tri 
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> for inlet with velocity ratio of 0.7 
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Jipira 13. - Continued. ‘Dilokoooa of ohadow aooe along innar inlet wall. Inlet reloolty ratio, 0.7) angle of attaok, 0°. 


CD 


MCA TH 3593 



4 .6 .8 1 2 

Ttaeijrocal of inertia 


Figure 14. - Outer-vall irajlngetaeut 










